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ABSTRACT 

The  foundation  of  the  theory  of  functionally  graded  plates  with  simply  supported  edges,  under  a 
Friedlander  explosive  air-blast,  are  developed  within  the  classical  plate  theory  (CPT).  Within  the  development 
of  the  theoiy,  the  two  constituent  phases,  ceramic  and  metal,  vary  across  the  wall  thickness  according  to  a 
prescribed  power  law.  The  theory  includes  the  geometrical  nonlinearities,  the  dynamic  effects,  compressive 
tensile  edge  loadings, and  the  damping  effects.  Also  presented  are  the  analytical  expressions  for  the  stresses,  in 
which  the  results  of  which  are  not  presented  here.  The  static  and  dynamic  solutions  are  developed  leveraging  the 
use  of  a  stress  potential  with  the  Extended-Galerkin  method  and  the  Runge-Kutta  method.  The  analysis  focuses 
on  how  to  alleviate  the  effects  of  large  deformations  through  proper  material  selection  and  the  proper  gradation 
of  the  constituent  phases  or  materials. 


1.  INTRODUCTION 

During  combat  situations,  the  structure  of  army  military 
vehicles  may  have  to  structurally  endure  the  effects  of 
blast  loading.  Advances  in  functionally  graded  materials 
(FGM’s)  which  combine  the  properties  of  two  dissimilar 
materials  has  been  a  motivating  factor  in  viewing  these 
types  of  materials  as  a  viable  alternative  to  the  current 
isotropic  metallic  structures  being  utilized  in  the  hull  and 
armor  plating.  FGM’s  are  microscopically 
inhomogeneous  with  thermo-mechanical  properties  which 
vary  smoothly  and  continuously  from  one  surface  to 
another.  These  graded  structures  allow  the  integration  of 
dissimilar  materials  like  ceramic  and  metals  that  combine 
different  or  even  incompatible  properties  such  as  hardness 
and  toughness. 

In  this  paper,  the  foundation  of  the  nonlinear  theory  of 
functionally  graded  plate-type  structures  under  an 
explosive  air-blast  is  developed.  An  approximate  solution 
methodology  for  the  intricate  nonlinear  boundary  value 
problem  is  devised,  and  results  that  are  likely  to 
contribute  to  a  better  understanding  of  the  structural 
behavior,  under  an  explosive  blast  with  beneficial 
implications  towards  their  improved  design  and 
exploitation  are  presented.  It  is  the  author’s  intent,  within 


this  paper,  to  fill  in  some  major  gaps  currently  existing 
within  the  specialized  literature. 

2.  BASIC  ASSUMPTIONS  AND  PRELIMINARIES 

The  plate  mid  surface  is  referred  to  a  cartesian 
orthogonal  system  of  coordinates  (x,y,z),  while  z  is  the 
thickness  coordinate  measured  positive  in  the  upwards 
direction  from  the  mid-surface  of  the  plate  with  h  being 
the  uniform  plate  thickness  and  y  is  directed  perpendicular 
to  the  x-axis  in  the  plane  of  the  plate.  See  figure  1  below. 


Figure  1:  A  simply  supported  functionally  graded  plate 
Shown  in  2-D,  under  an  explosive  blast. 
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The  nonlinear  elastic  theory  of  Functionally  Graded  (FG) 
Plates  is  developed  using  the  classical  plate  deformation 
Theory  [6],  It  is  also  assumed  that  the  FG  plate  is  made- 
up  of  ceramic  and  metal  phases  whose  material  properties 
vary  smoothly  and  continuously  across  the  wall  thickness. 
By  applying  the  rule  of  mixtures,  the  material  properties 
such  as  Young’s  Modulus,  Density,  and  Poisson’s  Ratio 
are  assumed  to  vary  across  the  wall  thickness  as 


Also,  for  the  purposes  of  simplicity  the  Poisson’s  ratio 
will  be  assumed  to  be  constant  throughout  the  plate 
structure. 

3.  KINEMATIC  EQUATIONS 
3.1  The  3-D  Displacement  Field 

Consistent  with  the  classical  plate  theory  [6],  the  distribution  of  the  3-1 


r(z)=rcVc(z)+rmvm(z),  (i) 

In  which  Pc  and  Pm  denote  the  material  properties  of  the 
ceramic  and  metallic  phases,  of  the  plate,  respectively. 

Vc(z)  and  Vm(z)  are  correspondingly,  the  volume 
fractions  of  the  ceramic  and  metal  respectively,  fulfilling 
the  relation 


Vc{z)+Vm(z)  =  1.  (2) 

By  virtue  of  (2),  Eq.  (la)  can  be  expressed  as 


/  \  /  \  dwn(x,  y,  t) 

u[x,y,  z,t)  =  u0(x,  y,t)-z — - , 

ox 

(  \  (  \  dWn\X,y,t) 

v{x,y,z,t)=v0{x,y,t)-z - - - ,  (6a-c) 

dy 

w{x,  y,z,t)=  w0(x,  y,  t) 


Within  these  equations,  ( u,v,w )  are  the  3-D 
displacement  quantities  along  the  (x,  y,  z)  directions, 
respectively.  While,  u0,  v0,  and  w0  are  the  2-D 
displacement  quantities  of  the  points  on  the  mid-surface. 


P(z)  =  (Pc-Pm)Vc(z)  +  Pm.  (3) 

Two  Scenarios  of  the  grading  of  the  two  basic 
component  phases,  ceramic  and  metal,  through  the  wall 
thickness  are  considered. 


3.2  NonLinear  Strain  Displacement 
Relationships 

The  nonlinear  strain  displacement  relationships  across 
the  plate  thickness  at  a  distance  from  the  mid-surface  take 
the  form 


Case  (1):  The  phases  vary  symmetrically  through  the 
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Where  the  signum  function,  sgn(z)=l,  z>0,  0,  z=0,  and  -1, 
z<0. 


N  is  termed  the  volume  fraction  index  which  provides  the 
material  variation  profile  through  the  plate  wall  thickness, 
(0  <  N  <  oo). 

Case  (2):  The  phases  vary  non-symmetrically  through 
the  wall  thickness,  and  in  this  case  there  is  full  ceramic  at 
the  outer  surface  of  the  plate  wall  and  full  metal  at  its 
inner  surface.  For  this  case,  Vc(z )  can  be  expressed  as 


Vc(z)  = 


(h  +  2z\ 


(5) 


It  should  be  noted  that  in  contrast  to  case  (2),  where 
there  exists  coupling  between  stretching  and  bending, 
such  coupling  is  not  present  for  the  symmetric  case  (1). 


Where, 
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(7) 


In  the  above  expressions,  ^£xx,£yy,Yxy),  are  referred  to  as 
the  membrane  strains  and  £yy>  Yxy )  are  referred  to  as 

the  flexural  bending  strains  which  are  also  known  as  the 
curvatures. 

4.  CONSTITUTIVE  EQUATIONS 

The  stress-strain  relationships  for  a  state  of  plane  stress 
is  expressed  as  [10] 
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The  material  stiffnesses,  Qy  (z ),  (i  =  1,2,6)  are  given  by 
[9,10] 
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The  standard  force  and  moment  resultants  of  a  plate  are 
defined  as 


N..,M..  \  = 
U  iJ 


f 


hi  2 


-A/2 


i(l,z)dz,  ( i,  j  =  x,y,xy  )  (10) 


With  the  use  of  Equations  (7)  -(10),  the  stress  resultants 
and  stress  couples  are  related  to  the  strains  by  [3] 


5.  GOVERNING  EQUATIONS 

Hamilton’s  principle  is  used  to  derive  the  equations  of 
motion  and  the  boundary  conditions.  It  is  formulated  as 

5J  =  5\\u+V-K)dt  =  Q,  (16) 

Where  t0  and  tqare  two  arbitrary  instants  of  time.  U 
denotes  the  strain  energy,  V  denotes  the  work  done  by 
surface  tractions,  edge  loads,  and  body  forces,  and  K 
denotes  the  kinetic  energy  of  the  3-D  body  of  the 
structure,  while  S  is  the  variational  operator.  In  Equation 
(16),  the  strain  energy  is  given  by 
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In  which  [A],  [B],  and  |  D  |  are  the  respective  in-surface, 
bending-stretching  coupling,  and  bending  stiffnesses.  For 
the  case  of  symmetric  FG  Plates,  [B]  =  0,  since  there  is 
no  bending-stretching  coupling.  The  global  stiffness 
quantities,  Ay,  By,  and  By,  (z,/=l  ,2,6)  are  defined  as 

Qvll,z,z2)dz,  (q./- 1,2,6)  (12) 

In  view  of  Equations  (3),  (4),  (5),  (9),  and  (12),  the 
global  stiffness  quantities  are  determined  as 

DA  1 ,  A2 ),  (A  1 ,  -®22 ),  (A  1 A22 )] =  J  {E 1  ’E2,E3), 

1-v 


(17) 


where  fl0  denotes  the  mid-surface  area  of  the  panel.  The 
work  done  by  external  loads  is  expressed  as 


SV  =  -  j"  Pt  (x,  y)Su^x,  y,^j-Ph  {x,  y)d^x,  y-  ^ 
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(18) 


In  the  above  expression,  Pt(x,  y)  is  the  distributed  force 
at  the  top  surface  (z  =  A/2),  Pb{x,y )  is  the  distributed 
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force  at  the  bottom  surface  (z  =  —h/ 2), 
iaxx’  ayy,axy  ayx)  are  the  specified  stress  components 
along  the  plate  edges,  and  ( Su ,  8 v)  are  the  virtual 
displacements  along  the  normal  and  tangential  directions, 
respectively,  along  the  plate  edges.  Considering  only  the 
transversal  inertia  of  the  structure,  the  kinetic  energy  is 
given  by 

K  =  \\p^)W2dV ,  (19) 


For  the  case  of  all  edges  simply  supported  and  freely 
movable  the  boundary  conditions  are  as  follows: 

=M„=  Nxy  =  0,  N„  =  N •„  on  x  =  0,  Lx 

w0  =  Myy  =  N„  =  0,  Nyy  =  Nyy  on  y  =  0,  L2  (22) 

It  should  be  mentioned  for  clarification  sake  that  for 

compressive  edge  loading  N*x  =  —Nxx,  and  N*y  = 

-N° 

nyy- 

6.  AIR-BLAST  LOADING 


Which  implies  that  the  variation  in  kinetic  energy  is 
expressed  as 


h 

SK=\  [2  p{z )WS\VdzdD.  0  (20) 

Jn0 

where  p(z)  is  the  mass  per  unit  volume. 

Considering  Equation  (16)  in  conjunction  with  Equations 
(17)-(20)  ,  along  with  the  constitutive  equations  (8),  the 
strain-displacement  relationships,  equations  (7),  and 
carrying  out  the  integration  through  the  thickness, 
integrating  by  parts  whenever  feasible,  using  the 
expression  of  global  stress  resultants  and  stress  couples, 
while  retaining  only  the  transversal  load,  transverse 
inertia,  and  transverse  damping,  invoking  the  arbitrary 
and  independent  character 

of  variations 

Su0,  8v0,  Sw0, 8{dw0/ dx)  ,  and  8{dw0/ dy)  one 

obtains  the  equations  of  motion  and  as  a  by-product  the 
boundary  terms  or  conditions.  This  results  in  three 
equations  of  motion  in  terms  of  the  global  stress  resultants 
and  stress  couples  and  the  four  boundary  conditions 
along  the  plate  edges.  These  equations  of  motion  and 
boundary  conditions  can  be  respectively  expressed  as 


With  the  ever  increasing  demands  for  increased  safety 
for  the  solders  in  the  field  to  operate  structurally  sound 
vehicles  in  the  event  of  an  improvised  explosive  device 
(IED)  or  some  other  type  of  explosive,  it  is  imperative 
that  an  understanding  of  the  structural  response  of  various 
components  within  military  combat  vehicles  under  an 
explosive  blast  be  understood  so  that  measures  can  be 
taken  from  a  design  standpoint  to  ensure  the  durability 
and  survivability  of  these  components.  To  begin  to 
achieve  this  understanding,  the  type  of  explosive  loading 
considered  here  is  a  free  in-air  spherical  air  burst.  Such  an 
explosion  creates  a  spherical  shock  wave  which  travels 
radially  outward  in  all  directions  with  diminishing 
velocity.  The  mathematical  expression  representing  the 
incident  blast  wave  from  a  spherical  charge  is  given  by 


Pt(t)={Pso-Po)  1 


t-tn 
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Pso  is  the  peak  overpressure  above  ambient  pressure,  P0  is 
the  ambient  pressure,  ta  is  the  time  of  arrival,  tp  is  the 
(21a)  positive  phase  duration  of  the  blast  wave,  and  t  is  the 
time.  The  waveform  shown  in  figure  4  is 
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In  equations  (23)  and  (24)  Z  is  known  as  the  scaled 
distance  given  by  Z  —  R/W1^3  with  R  being  the  standoff 
distance  in  meters  and  W  being  the  equivalent  charge 
weight  of  TNT  in  terms  of  kilograms.  Also,  a  is  known  as 
the  decay  parameter  which  is  determined  by  adjustment  to 
a  pressure  curve  from  a  blast  test. 

For  the  conditions  of  standard  temperature  and  pressure 
(STP)  at  sea  level,  the  time  of  arrival  ta,  and  the  positive 
phase  duration  tp,  can  be  determined  from  [4] 
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-(1-U)V2AT 


(29) 


Where  t1  represents  either  the  arrival  time  or  positive 
phase  duration  for  a  reference  explosion  of  charge  weight 
Wl7  and  t  represents  either  the  arrival  time  or  positive 
phase  duration  for  any  explosion  of  charge  weight  W .  The 
determination  of  the  standoff  distance  for  any  charge 
weight  W  follows  a  similar  reasoning.  The  application  of 
these  relationships  is  known  as  cube  root  scaling.  It 
should  be  understood  that  in  applying  these  relationships 
that  the  standoff  distances  are  themselves  scaled 
according  to  the  cube  root  law. 

7.  SOLUTION  METHODOLOGY 

To  satisfy  the  first  two  equations  of  motion,  equations 
(31a,b),  a  stress  potential  will  be  utilized  which  allows  the 
in-plane  stress  resultants  to  be  expressed  by  letting 


N  =  d  j\f 
xx  2  ’  y 
dy 


dl±  N 

fix2  ’  **  dxdy 


(26) 


The  third  equation  of  motion,  equation  (21c),  can  be 
expressed  in  terms  of  two  unknown  variables,  the  stress 
potential  0  and  the  transverse  displacement  w0  which 
results  in 


DV  w0  - 
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Where, 
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In  equations  (27)  and  (29),  V2=  d2/dx2  +  d2/dy2 

where  V  is  referred  to  as  the  Laplacian  operator. 

Equations  (27a)  and  (29)  are  the  basic  governing 
equations  used  to  investigate  the  structural  response  of  FG 
plates  under  external  excitation  loading.  For  the  purposes 
of  this  paper,  from  this  point  forward,  the  thermal  terms 
will  be  discarded.  To  this  end,  to  solve  equations  (27a) 
and  (29),  the  approach  adopted  from  [2]  will  be  utilized. 
In  this  respect,  the  following  functional  forms  are 
assumed  for  wc  and  <p  [2]. 


w0(x,y,t)  =  wmn  (t)  sinAmx  sin /jny 

(30a, b) 

<t>  =  Amn  (')cos  2Amx  +  (i)cos  2  nny  + 

Cmn  (t )  cos  2Xmx  cos  2  iuny  +  Dmn  (t)sin22mxsin2//„y  + 


{Key2  +\N;yx2 


Where  Am  —  rrm/Lx,\in  —  nn/L2:  m,n  —  1,2, ...  are  the 
number  of  half  waves  in  the  x  and  y  directions, 
respectively,  and  wmn(t)  is  the  amplitude  of  deflection. 
Also,  d-jun  (l),  drnn  (t) ,  Cmn  (t),  and  Dmn  (t)  are 

coefficients  to  be  determined.  By  substituting  equations 
(50a, b)  into  equation 

(49),  the  coefficients 

Cmn (t),  and  Drnn(t )  are  determined  as 


Anm(t)  = 


E\  w~„ 


MV2,  _  /\  Exw2mn{f)y l2 


32A; 


"  •>  B m„  (t ) 


32  Hi 


(31) 


Cmn  (t)=Dnm(t)=Q 


By  letting, 


This  gives  one  governing  equation  with  two  unknowns, 
w0  and  0.  One  more  equation  is  needed  in  terms  of  the 
same  unknowns  which  will  give  two  equations  in  terms  of 
two  unknowns  which  can  then  be  solved.  This  will  come 
from  the  compatibility  equation.  By  eliminating  the  in¬ 
plane  displacements  from  the  strain-displacement 
relationships,  equations  (7)  the  relationship  between  the 
in-plane  strains  and  the  transversal  deflection  known  as 
the  compatibility  equation  can  be  shown  to  be  given  by 


g24 

dy2 


d2s° 
u  byy 


d2r° 


xy 


dx2  dxdy 


f  p2  \2 

o  w0 
y dxdy ^ 


5° w0  d2w0 
dx2  d2y 


(28) 


Emn  (t)sinsin  A,mxsin /Jny  (32) 

And  integratating  both  sides  over  the  plate  area  gives 


«)  =  ^>,  (m,n)  =  (1, 1)  (33) 

71 

and  introduction  of  equations  (30a, b)  and  (31)  into 
equation  (27a)  and  retaining  the  resulting  equation  along 
with  the  unsatisfied  boundary  conditions  in  the  energy 
functional  and  applying  the  Extended  Galerkin  technique 
results  in  the  following  nonlinear  differential  equation 


Which,  in  terms  of  0  and  w0,  becomes 
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governing  the  structural  response  of  FG  plates,  under 
external  excitation. 


Nxx  ~A\sxx  £yy  +  \sxx  +  B12 £yy 


W mn  W  ^^mn^mnAmn  W  +  ^mn^mn  W  mn  W 

—  Bmn  (0 

(34) 

Where,  wmn(t)  represents  the  amplitude  of  deflection  of 
the  plate  as  a  function  of  time  The  expressions  for 

<umn,  Amn,  flmn,  ancL  Pmn  are  not  provided  here. 

Equation  (34)  is  a  nonlinear  equation  in  terms  of  the 
plate  deflections  as  a  function  of  time.  To  obtain  the  plate 
deflections  as  a  function  of  time,  equation  (34)  is  solved 
using  the  Fourth-Order  Runge-Kutta  Method. 

8.  MECHANICAL  STRESS  ANALYSIS 

With  the  transversal  deflections  in  hand,  as  a  function  of 
time,  the  plane  stresses  can  be  determined  from  the 
constitutive  equations  (8)  with  the  use  of  the  strain- 
displacement  relationships,  equations  (7)  and  the  stiffness 
equations,  equation  (9)  results  in. 


E{z) 


® XX  "  9 

1-v2 


du, 

dx 


o  +VA  +  I(  Swo 


dy  2\  dx 


v 

+  — 
2 


cWq 

~dy 


g-WQ  |  v  d-WQ 


dxz 


gy 


(35a) 


Gyy  ~ 


E(z) 


1  —  v 


dvn  dll(]  1 

— -  +  v — -  +  — 
dy  dx  2 


dwQ 

~dy 


+  vf  dwt 0 


2  1  dx 


d  w0  |  ^  d  w0 


axy  ~ 


dyz 

E(z) 


dxz 


2d  +  v) 


duQ  Dvq  Dwq  dwn  DzWq 
4 - 1 - 2  z  • 


dy  dx  dx  dy 


dxdy 


(35b) 


(35c) 


In  order  to  determine  the  stresses,  the  expressions  for 
no  and  Vo  in  terms  of  the  transversal  amplitude  of 
deflection,  wmn  ( t )  needs  to  be  determined.  These 
relationships  can  be  determined  from  the  global 
constitutive  equations  (11).  Expanding  equation  (11) 
gives, 


Nyy  -  A\2£xx  +  A22£yy  +  B\2£xx  +  ^22 £yy  (36a, b) 


Solving  Equations  (36a, b)  simultaneously  for 

,  and  in  terms  of  Nxx,  Nvv,  £xx,  and  s ^  and 
utilizing  equations  (7)  and  (26)  gives, 


0vo_ 

gy 


/  \ 

A\ 

d2<P  + 

(  \ 

-An 

d2<p 

2 

VA12  -AnA2 2  ) 

dx2 

2 

K^12  ~A11A22  ) 

dy 2 

AnBn  ~  Ai  i  g|2 

|  g  2  w0  | 

/  \ 

^12^12  ~  ^11^22 

^2 

O  Wq 

2 

V  A12~A11A22  ) 

<N 

c§ 

2 

Hi2-HiiH22  ) 

dy 2 

''dw, o 

dy  j 


1 

o  I 
CO  | 

(  \ 

~A\2 

d2(p 

_ _  _|_ 

( 

A22 

d2<p 

dx 

k,A12  ~AUA22  j 

dx2 

\A12 

~  A\  \A22  j 

dy2 

f  \ 

A\2B\2  ~A22B\l 

^2 

O  Wq 

/  \ 

^12^22  ~A22B12 

^2 

O  Wq 

v  An~AUA22  j 

dx2 

v  ^12  —  ^1 1^22  j 

dy2 

f  dw,  A 


dx 


(37a, b) 


With  the  use  of  equations  (13a-c),  (30a, b),  and  (31)  and 
integrating  equations  (37a, b),  z<0and  v0  can  be 
determined.  Having  in  hand,  the  expressions  for  u  () ,  v0  , 
and  w’o  ,  the  expressions  for  the  stress  distribution  can  be 
determined  from  equations  (35a-c). 


9.  RESULTS  AND  DISCUSSION 


To  illustrate  the  present  approach,  a  ceramic-metal 
functionally  graded  plate  consisting  of  Ti-6A1-4V  and 
Aluminum  Oxide  with  the  following  material  properties, 
which  were  adopted  from  [9],  were  considered  for  the 
numerical  results  presented. 


Ec  =  320.24  GPa, 
Em  =105.7  GPa, 


pc  =3750  kg/m3, 
Pm  =4429  kg  /m3, 


vc  =  0.26 
vm  =0.2981 


The  geometrical  properties  used  for  the  FG  Plate  are 
Lx  =  lm,  (//  =  Z,1/Z,2  =1.25,  and  unless  otherwise  stated 
h  =  0.0254  m  and  the  halfwaves,  (m,ri)  =  (1,1).  In 
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addition,  the  following  reference  values,  in  Table  1.  were 
utilized  to  determine  the  time  of  arrival  and  positive  phase 
duration  [5]. 


Table  1.  Airblast  Parameters  Versus  Distance 
for  a  One  Kilogram  (lVt )  TNT  Spherical  Air 
Burst  [5]. _ 


Standoff 
Distance,  R1 
(m) 

Arrival 
Time,  tai 
(msec) 

Positive  Phase 
Duration,  tpi 
(msec) 

1.0 

0.532 

1.79 

In  Fig.  2  below,  comparisons  are  made  between  the 
deflections  of  an  asymmetric  and  symmetric  plate.  It  can 
be  seen  that  the  deflections  are  a  little  out  of  phase  with 
each  other.  Also,  the  symmetric  FG  plate  has  slightly 
lower  deflections  for  a  volume  fraction  index  of  0.1.  As 
can  also  be  seen,  that  for  a  fixed  amount  of  damping,  the 
deflections  are  attenuated  rather  quickly. 


Fig.  2  .  Implications  of  the  symmetiy  on  the  central  deflection  as  function  of 
time.  (A/=0.1,  Aj  i=0.1) 


in  Fig.  3,  It  is  shown  for  an  asymmetric  FG  plate,  that  as 
the  volume  fraction  index  increases  from  a  value  of  zero 
(fully  ceramic)  to  a  value  of  10  (essentially  fully  metal) 
that  the  deflections  increase  then  taper  off  as  phases  of  the 
constituent  materials  approach  fully  metal.  This  shows 
that  metal  have  higher  deflections  than  the  ceramic 
implying  that  the  grading  of  the  FG  plate  plays  an 
important  role. 

In  Fig.  4,  which  is  the  counterpart  of  Fig.  3  for  a 
symmetric  FG  plate,  reveals  similar  behavior.  As 
mentioned  earlier  symmetric  FG  plates  have  lower 


deflections  than  asymmetric  FG  plates.  Also,  it  can  be 
seen  that  damping  attenuates  the  deflections  very  rapidly. 


Fig.  3  .  Implications  of  the  volume  fraction  index  on  the  central  dflcction  as  a  function  of 
time  on  an  asymmetric  functionally  graded  plate.  (A|  j=0. 1) 


Fig.  4.  Implications  of  the  volume  fraction  index  on  the  central  dflcction  as  a  function  of 
time  on  an  symmetric  functionally  graded  plate.  (A|  |=0. 1) 

10.  CONCLUDING  REMARKS 

A  rigorous  treatment  of  functionally  graded  plates  with 
grading  in  the  transverse  direction  has  been  studied.  It  has 
been  shown  that  damping  has  an  important  effect  when  it 
comes  to  the  attenuation  of  the  deflections.  It  has  also 
been  shown  that  the  symmetry  of  the  transverse  grading 
throughout  the  structure  plays  an  important  role  in  the 
deflection-time  history  of  the  structure. 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 
Page  7  of  8 


From  a  design  standpoint,  it  would  be  appropriate  at  this 
point  to  state  that  integration  of  functionally  graded 
materials  within  plate-type  structures  would  benefit  the 
structural  response  of  the  structure.  Also,  it  should  be 
mentioned,  although  not  shown  here,  that  the  choice  of 
the  ceramic  and  metal  constituent  materials  chosen  would 
also  have  a  great  impact  on  the  response  of  the  structure. 

The  idea  is  to  reduce  the  stresses  within  the  structure 
concerned  here.  By  reducing  the  magnitude  of  the 
deflections,  the  stresses  are  reduced.  It  is  hoped  and 
realized  that  this  present  study  presented  here  will  give 
insight  into  some  of  the  factors  that  can  play  an  important 
role  in  the  structural  response  of  functionally  graded 
plates  and  fill  in  some  of  the  fundamental  missing  gaps 
within  this  subject  area. 

FUTURE  WORK 

It  should  be  noted  that  further  work  should  and  needs  to 
be  explored  which  would  address  comparing  finite 
element  blast  modeling  and  simulation  results  with  the 
current  analytical  results  based  on  the  theory  of  elasticity. 
This  can  be  accomplished,  by  exploring  the 
implementation  of  these  analytical  equations,  as  a  user 
subroutine,  in  one  of  the  applicable  software  tools. 
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